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ABSTRACT 

We present new results based on observations carried out with the Rossi X-ray Timing 
Explorer during the decay of an outburst of the low-mass X-ray binary (LMXB) and atoll 
source 4U 1608-52. Our results appear to resolve, at least in 4U 1608-52, one of the long- 
standing issues about the phenomenology of the kilohertz quasi-periodic oscillations (kHz 
QPOs), namely, the lack of a unique relation between the frequency of the kHz QPOs 
and the X-ray flux. We show that despite its complex dependence on the X-ray flux, 
the frequency of the kHz QPOs is monotonically related to the position of the source 
in the color-color diagram. Our findings strengthen the idea that, as in the case of Z 
sources, in the atoll sources the X-ray flux is not a good indicator of M, and that the 
observed changes in the frequency of the kHz QPOs in LMXBs are driven by changes in 
M. These results raise some concern about the recently reported detection of the orbital 
frequency at the innermost stable orbit in 4U 1820-30. 
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1. Introduction 

Nearly three years have elapsed since the first kilo- 
hertz quasi-periodic oscillations (kHz QPOs) were 
discovered with the Rossi X-ray Timing Explorer 
(RXTE) i n the X-ray flux of Scorpius X-l (van der 
Klis et al. jL996|) and 4U 1728-34 (Strohmayer, Zhang, 
& Swank |1996| ). In the meantime, kHz QPOs have 
been observed in the persistent flux of 1 8 low mass 
X-ray binaries (LMXBs; see van der Klis 1998 for a 
review), both in the so called Z sour ces a nd in the 
atoll sources (Hasinger & van der Klis 1989 , hereafter 
HK89). Except in Aql X-l, which showed a single 
kHz peak in its power spectrum (Zhang et al. 1998a), 



in the other 17 sources two simultaneous kHz QPO 
peaks have been observed. In some sources a third, 
nearly-coherent, QPO peak has been detected dur- 
ing type-I X-ray bursts, with a frequency which was 
consistent with being equal to one or two times the 
frequency separation of the kHz QPOs observed in the 



persistent flux (Strohmayer, Swank, & Zhang 1998). 
It was suggested that a beat frequency mechanism 
is responsible for this commensurability in the QPO 
frequencies (Strohmayer et al. 1996; Miller, Lamb & 
Psaltis 1998), but this interpretation is not without 
proble ms (van der Klis et al. 1997; Mend ez et a l. 
1199801 ; Mendez, van der Klis, & van Paradijs |l998b| ). 

The dependence of the kHz QPO frequencies on 
X-ray luminosity, which is usually assumed to be a 
measure of the mass accretion rate, is complex. While 
in a given source, on a time scale of hours, there is 
a good correlation between frequency and luminosity, 
sources that span nearly three orders of magnitude in 
luminosity, such as ScoX-1 and 4U 0614+09, show 
kHz QPOs th at co ver the same range of frequencies 
(van der Klis 1997 ). It is as if the frequency of the 
kHz QPO depends on the difference between instan- 
taneous and average luminosity in each source rather 
than on the luminosity itself. 

A similar effect is seen between observations of the 
same source at different epochs. On time scales of 
hours, frequency and X-ray flux are well correlated, 
but between different epochs the source covers the 
same range of frequencies even if the average flux is 
differe nt by 40 % or more (e.g., Aql X-l; Zhang et al. 
1998a| ). 

In this paper we present new results that appear 
to resolve the latter of these two issues, at least in 
4U 1608-52. We show that while on time scales longer 
than ~ 1 day the frequency of the kHz QPOs is not 



well correlated to the X-ray flux, it is very well corre- 
lated to the position of the source in the color-color 
diagram. From this result we conclude that the ob- 
served changes in the frequency of the kHz QPOs in 
4U 1608-52 are driven by changes in the mass accre- 
tion rate, and that the lack of correlation between 
QPO frequency and X-ray count rate occurs because, 
as in the case of the Z sources, in atoll sources there 
is no one-to-one relationship between the observed X- 
ray flux and the mass accretion rate. 

2. Observations and Results 

All the observations presented here were obtained 
using the proportional counter array (PCA) onboard 
RXTE. We include the data o f the decay of t he 199 6 
outburst (Berger et al. 1996; Mendez et al. 1998a), 



and of the 1998 outburst (Mendez et al. |l998c| ; Fig. 1 
there shows a light curve of the 1998 outburst.) The 
observations as well as the modes used to record the 
data are described in Mendez et al. (1998c). We also 
include here a recent Public Target of Opportunity 
RXTE/PCA observation of ~ 7.3 ks performed on 
August 6 1998, 03:13 UTC. The observing modes for 
this last observation were similar to those used by 



Mendez et al. (1998c) after 1998 March 27. 



We calculated count rates in 5 energy bands, 2.0 — 
3.5 - 6.4 - 9.7 - 16.0 keV, and 2.0 - 16.0 keV, tak- 
ing into account the gain changes applied to the PCA 
in March and April 1996. In a few of the observa- 
tions one or two of the five detectors of the PCA were 
switched off; we only used the three detectors which 
were always on to calculate these count rates. We 
subtracted the background contribution in each band 
using the standard PCA background model version 
2.0oP], and normalized the count rates to 5 detectors. 

In Figure [l] we show a color-color diagram of 
4U 1608-52. This is the first time that 4U 1608-52 
is observed to move across all the branches of the 
atoll, and constitutes one of the best examples of the 
color-color diagram of an atoll source. Based on this 
diagram we conclude that in 1996, as the source count 
rate decreased, 4U 1608-52 gradually moved from the 
lower part of the banana to the island state. In 1998, 
at the peak of the outburst, 4U 1608-52 was in the 
upper part of the banana, and gradually moved down 
to the lower part of the banana and the island state as 



The PCA Background Estimator is available at 



littp: / /lheawww. gsfc. nasa.gov / users /stark/pea/ pcabackest . html , 
which is maintained by N ASA/GSFCJ. 



2 



the count rate decreased. In general, the count rate is 
observed to increase along the track from the island 
state to the banana branch, but the relation between 
count rate and either of the two colors is much less 
clean than that between colors. 

The high-time resolution data confirm this prelim- 
inary state classification. We divided the 2 — 60 keV 
data into segments of 256 s and 512 s, and calculated 
power spectra up to a Nyquist frequency of 2048 Hz, 
normalized to fractional rms squared per Hertz. The 
characteristics of the <100 Hz part of these power 
spectra changed in correlation with the position of 
the source in the color-color diagram. When 4U 1608- 
52 was in the upper and lower parts of the banana, 
the power spectra fitted a power law below ~ 1 Hz 
(the Very Low Frequency Noise, VLFN), and an ex- 
ponentially cut-off power law above ~ 1 Hz (the High 
Frequency Noise, HFN; see HK89). As 4U 1608-52 
moved from the upper to the lower parts of the ba- 
nana, the fractional amplitude of the VLFN (0.001 — 1 
Hz) decreased from ~ 6 % rms to ~ 2 % rms, while 
the fractional amplitude of the HFN (1 — 100 Hz) in- 
creased from ~ 1 % rms to ~ 5 % rms. In the island 
state the VLFN disappeared completely (the 95% 
confidence upper limits were ~ 1 % rms), and the am- 
plitude of the HFN increased further to ~ 10 - 17% 
rms. (A more detailed analysis of the low frequency 
part of the power spectra will be presented elsewhere.) 

The >100 Hz part of the power spectra also changed 
in correlation with the position of the source in the 
color-color diagram: we only observed kHz QPOs 
when 4U 1608-52 was at certain positions of the color- 
color diagram (see Fig. |l|). For those segments where 
we observed QPOs, the 2 — 60 keV fractional am- 
plitudes of the lower-frequency and higher-frequency, 
hereafter the lower and upper QPO, varied from 5.3 % 
to 9.1% rms, and from 3.3% to 8.8% rms, respec- 
tively. We did not detect kHz QPOs in the upper 
part of the banana, with 95 % confidence upper limits 
of 0.8 % to 4.6 % rms depending on the source count 
rate and the assumed width of the QPO, nor in the 
extreme island state (upper right corner of the color- 
color diagram), with 95% confidence upper limits of 
3.5% to 10% rms. While these upper limits strongly 
suggest the absence in the upper banana of kHz QPOs 
as strong as those observed in the lower banana and 
the island states, we cannot rule out the presence of 
such kHz QPOs in the extreme island state, when the 
count rates were lowest. 

To further characterize the dependence of the kHz 



QPOs on other source parameters, we selected only 
those data where we detected two simultaneous kHz 
QPO peaks in the power spectrum (see Mendez et 
al. 1998c), so that the identification of the observed 
kHz peaks is unambiguous. We divided the data in 
segments of 64 s and produced a power spectrum for 
each segment extending from 1/64 Hz to 2048 Hz. In 
Figure ^| we show the dependence of the frequency 
of the lower QPO, ^low, a s a function of count rate 
for 4U 1608-52. This figure shows several branches, 
which reflect the relation between v\ ow and count rate 
during individual observations that span from ~ half 
an hour to ~ 8 hours. The only exception is the 
branch at the lowest count rate where two different 
observations taken 8 days apart overlap. This figure 
clearly shows that while on time scales of a few hours 
or less flow is well correlated to count rate, on time 
scales greater than ~ l day the relation is complex 
and flow is not uniquely determined by the count rate. 
We obtained the same result using the 2 — 60 keV 
source flux instead of the count rate. 

There is a much better correlation between v\ ovl 
and the position of the source on the color-color di- 
agram. In Figure || we show v\ ow as a function of 
hard color (see Fig. |l|) for the same intervals shown 
in Figure |2| The complexity seen in the frequency 
vs. count rate diagram (Fig. ||) is reduced to a single 
track in the frequency vs. hard color diagram. The 
frequency of the lower kHz QPO increases as the hard 
color decreases, i.e., as 4U 1608-52 moves from the is- 
land state to the lower banana, and keeps increasing 
at the turn of the lower banana, where the hard color 
reaches its lowest value of ~ 0.40. 

The shape of the track in Figure |^ suggests that 
the hard color may not be sensitive to changes of state 
when the the source moves into the banana in the 
color-color diagram. To further investigate this, we 
applied the Sz parameterization (e.g., Wijnands et al. 



1997b), which we call S a for atoll sources, to the color- 



color diagram. In this phcnomcnological approach we 
approximated the shape of the color-color diagram 
with a spline (Fig. ^), and we used the parameter S a 
to measure positions along this spline. We set S a to 
1 at (2.67,0.77), and to 2 at (2.19,0.42), as indicated 
in Figure ^. We measured the position on the color- 
color diagram of each of the 64-s segments. We then 
grouped these segments into 36 sets, such that within 
each set v\ ow did not vary by more than ~ 25 Hz, 
and the source count rate did not vary by more than 
10%. (This is to avoid mixing segments that come 
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from different branches in Figure g.) For each set 
we calculated the average frequency of the lower and 
upper kHz QPO peak, and of S a . In Figure [|we plot 
the frequencies of both kHz QPOs as a function of 
S a . The error bars are the standard deviation of each 
selection. We also include in this figure several extra 
sets (open squares) for which we only detected one of 
the kHz QPOs, which therefore we could not a priori 
identify as the upper or lower peak. As expected, 
S a is more sensitive than the hard color to changes 
of state when the the source moves from the island 
state into the turn of the banana: S a increases from 
~ 2.04 to ~ 2.15 there, while the hard color saturates 
at ~ 0.40. Figure ^ confirms that the frequency of the 
kHz QPOs is strongly correlated to the position of the 
source in the color-color diagram. It also shows that 
the steep rise of kHz QPO frequency at the turn of the 
lower banana is not just due to a lack of sensitivity 
of the hard color to state changes in this part of the 
color-color diagram, but that it occurs as a function 
of the position in the color-color diagram as well. 

3. Discussion 



Previous observations of kHz QPOs in atoll sources 
have resulted in a confusing picture concerning the de- 
pendence of QPO frequency on mass accretion rate. 
Cases have been reported where, in a given source, 
QPO frequency showed a good correlation to count 
rate or spect ral hardness (For d, van der Klis, & 
Kaa ret 1998a ; Wijnan ds et a l 
al. 



1998| ; Ford et al. 1998b] Z hang et al. 1998t 



1997c; Wijnands et 



Wijnands & van der Klis 1997). In other cases, a 
correlation with count rate or flux was conspicuously 
lackin g (For d et al. |1997| ; Mendez et al. |l998e| ; Zhang 
et al. Il998a|). 



Our observations show for the first time that a total 
lack of correlation between frequency and count rate 
on time scales longer than a day (Fig. ||) can coexist 
with a very good correlation between frequency and 
position in the X-ray color-color diagram (Fig. ||). 
The frequency increases with S a , as the source moves 
from the island to the banana. Only on time scales of 
hours does the QPO frequency appear to also corre- 
late well with count rate. The presence of the QPOs 
also correlates well with the position in the color-color 
diagram: the QPOs are only detected in the lower ba- 
nana and the moderate island states, and disappear 
both in the upper banana and in the extreme island 
states (Fig. |). 



In atoll sources M is thought to increase monoton- 
ically with 5 a along the track in the color-color dia- 
gram, from the island to the upper banana (HK89), 
whereas X-ray count rate tracks M much less well 
(van der Klis et al 
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van der Klis 1994; Prins 
& van der Klis [19971 ). The properties of the < 100 
Hz power spectra depend monotonically on inferred 
M (HK89). Our result that the frequency of the 
kHz QPO is well correlated to S a , but not to X- 
ray count rate, implies that in 4U 1608-52 the kHz 
QPO frequency also depends monotonically on in- 
ferred M. By extension this conclusion applies to 
each atoll source; in Z sources similar conclusi ons wer e 
previously also reached (e.g., Wijnands et al. 1997a ). 
Our analysis, however, sheds no light on the question 
why sources with very different inferred M (e.g., 4U 
0614+09 and Sco X-l) can have kHz QPO in the same 
range of frequencies. 

Further indications for this interpretation come 
from the simultaneous analysis of the low and high 
frequency parts of the power spectra of other atoll 
sources. In 4U 1728-34 the kHz QPO frequencies 
were recently found to be very well correlated to sev- 
eral < 100 Hz power spectral properties (Ford & van 
der Klis 1998), and Psaltis, Belloni and van der Klis 
(1998) obtained a similar result for a number of other 
atoll (and Z) sources. In all these sources, not only 
the position in the color-color diagram and the vari- 
ous low frequency power spectral parameters, but also 
the frequencies of the kHz QPOs are all well corre- 
lated with each other. This indicates that the single 
parameter, inferred to be M, which governs all the 
former properties also governs the frequency of the 
kHz QPO. 

X-ray intensity is the exception: it can vary more 
or less independently from the other parameters. In 
4U 1608-52, it can change by a factor of ~ 4 (see 
Fig. |^) while the other parameters do not vary sig- 
nificantly. If as inferred, this constancy of the other 
parameters means that M does not change, then this 
indicates that strongly variable beaming of the X-ray 
flux, or large-scale redistribution of some of the ra- 
diation over unobserved energy ranges is occurring in 
order to change the flux by the observed factors, with- 
out any appreciable changes in the X-ray spectrum. 
We may need to scrutinize more closely the concept 
of M in order to solve this dilemma. For example, 
perhaps the M governing all the other parameters is 
the M through the inner accretion disk, whereas mat- 
ter also flows onto the neutron star in a more radial 
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inflow, or away from it in a jet. 

In 4U 1608-52 we observe no evidence for a satura- 
tion of the frequency of the kHz QPOs at a constant 
maximum value as M increases, different from what 
Zhang et al. fll998b|) inferred for 4U 1820-30. They 



presented data in which the kHz QPO frequencies in- 
crease with count rate up to a threshold level, above 
which the frequencies remain approximately constant 
while the count rate keeps increasing. Interpreting 
count rate as a measure for M they argue that this 
is evidence for the inner edge of the disk reaching 
the general-relativistic innermost stable orbit. How- 
ever, we have shown here that, at least in 4U 1608-52, 
count rate is not a good measure for M. Inspection of 
Figure ^| suggests that with sparser sampling our plot 
could easily have looke d similar to the one presented 
by Zhang et al. ( |l998b| ) for 4U 1820-30. It will there- 
fore be of great interest to see if in 4U 1820-30 the 
saturation of QPO frequency as a function of count 
rate is still there when this parameter is plotted as a 
function of position in the X-ray color-color diagram. 
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Fig. 1. — Color-Color diagram of 4U 1608-52. The soft and hard colors are defined as the ratio of count rates in the 
bands 3.5 - 6.4 keV and 2.0 - 3.5 kcV, and 9.7 - 16.0 keV and 6.4 - 9.7 kcV, respectively. The contribution of the 
background has been subtracted, but no dead-time correction was applied to the data (in this case the dead-time 
effects on the colors are less than 1 %). Each point represents 128 s of data. We show the typical error bars in the 
banana and the island states. Filled and open symbols indicate segments with and without kHz QPOs, respectively. 
We represent the island and banana states (defined on the basis of the low-frequency part of the power spectra) by 
circles and triangles, respectively. The curve shows the parametrization of the color-color diagram in terms of 5 a 
(see text). 
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Fig. 2. — Relation between the frequency of the lower kHz QPO and the 2 — 16 keV count rate. The count rates 
have been corrected for background, and normalized to 5 detectors. Each point represents a 64-s segment (see 
text). We only include data where both kHz QPOs are detected simultaneously , so that we can unambiguously 
identify the lower kHz peak. 
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Fig. 3. — Relation between the frequency of the lower kHz QPO and the hard color (see Fig. 1), for the same 
segments shown in Figure 2. 
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Fig. 4. — Diagram of the frequencies of the upper and lower kHz QPOs vs. the position of the source on the 
color-color diagram, as measured by 5 a (see Fig. 1, and text). Circles and filled squares represent the lower and 
the upper kHz peak, respectively. Open squares represent segments where we only detected one of the kHz QPOs, 
and therefore we could not identify it as the upper or lower peak. The asterisk is one observation where we see a 
single QPO at a level of significance of 3a only. 
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